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We present a new limit on the cosmological constant based on the absence of correlations between 
the cosmic microwave background (CMB) and the distribution of distant radio sources. In the 
cosmological constant-cold dark matter (ACDM) models currently favored, such correlations should 
have been produced via the integrated Sachs- Wolfe effect, assuming that radio sources trace the local 
{z 1) matter density. We find no evidence of correlations between the COBE 53Hz microwave 
map and the NVSS 1.4 GHz radio survey, and obtain an upper limit for the normalized cross- 
correlation of {SN5T)/y/{5N^}{ST2) < 0.067 at the 95% CL. This corresponds to an upper limit on 
the cosmological constant of Qa < 0.74, which is in marginal agreement with the values suggested by 
recent measurements of the CMB anisotropy and type-IA supernovae observations, Ha — 0.6 — 0.7. 
If the cosmological model does lie in this range, then the integrated Sachs- Wolfe eff'ect should be 
detectable with upcoming CMB maps and radio surveys. 
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Recent observations of supernovae light curves [y sug- 
gest that the expansion of the universe is accelerating 
rather than decelerating. Combined with evidence from 
the cosmic microwave background (CMB) and a 

number of other observations this suggests the uni- 
verse is spatially flat and dominated by a cosmological 
constant, SIa = A/SiJg ~ 0.6 — 0.7, where Hq is the 
Hubble expansion parameter. Such a low value of A is 
difficult to explain from fundamental grounds, so it is 
vital that we try to confirm this result by other means. 

CMB anisotropics can arise via the integrated Sachs- 
Wolfe effect (ISW) as the photons travel through 
the time-dependent gravitational potentials of collapsing 
structures. One consequence of a large cosmological con- 
stant is that such time-dependent potentials exist even 
on very large scales where the collapse is linear, which 
is not the case for a flat, matter dominated universe. 
These fluctuations are likely to be small compared to 
those imprinted at the surface of last scattering (redshifts 
z ^ 1000) and arc difflcult to detect directly, however, 
they can be observed by looking for spatial correlations 
between the CMB and the nearby matter density 0,1). 
This requires a probe of the matter density out to red- 
shifts of 2; 2 and suggested candidates include radio 
galaxies, quasars and the x-ray background. 

Other processes can also lead to correlations between 
the CMB and the local matter density. These include 
gravitational lensing, scattering from hot electrons (the 
Sunyaev-Zeldovich effect) and photons passing through 
the time-dependent potentials of non-linear collapsing 
structures (the Rees-Sciama effect). While the study of 
these effects can also beneflt from cross correlation analy- 
ses Q , the ISW effect is unique in that it occurs on very 
large scales {9 > 1°) where the fluctuations are simple 
and linear. 



In the first attempt to detect this effect, Boughn, Crit- 
tenden, and Turok |l^ cross-correlated the CMB with 
the hard (>2 keV) x-ray background, which is thought 
to arise primarily from active galactic nuclei out to a red- 
shift of z - 4 Jul . A COBE-DMR CMB map ^ was 
cross-correlated with the HEAOl x-ray map |l^ and no 
significant correlation between the two maps was found. 
The interpretation of this result is difflcult since there 
has been no unambiguous measurement of correlations 
in the x-ray background, so the x-ray bias factor, the ex- 
tent to which x-rays trace matter, is largely unknown. 
For a large x-ray bias, 4, the implication is that 

Ha < 0.5, in confiict with the currently favored A cold 
dark matter (ACDM) cosmological model; for no bias- 
ing, bx = I, the limit is much weaker: CIa < 0.95 [ p^ . 
The true bias may well be time dependent and will not 
be known until the structure in the x-ray background is 
definitively detected. 

Here we attempt to detect the ISW effect by cross- 
correlating the CMB with a deep radio source survey. 
While the Poisson noise due to the finite number den- 
sity of radio sources is relatively large, this analysis has 
the advantage that the clustering properties of the ra- 
dio sources have been measured and so the bias factor 
for these sources can be deduced within the context of a 
particular cosmological model. 

The small correlations introduced by the ISW effect 
can be contaminated by accidental correlations of the lo- 
cations of radio sources with the relatively large CMB 
fluctuations that originate at higher redshift. To mini- 
mize these accidental correlations, it is necessary to av- 
erage over many statistically independent regions of the 
sky. Since the ISW correlations are on angular scales 
of several degrees, this means using surveys that cover 
as much of the sky as possible, even if they have more 
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noise than smaller area surveys. Thus, in the radio we 
use the NVSS radio source survey (82% sky coverage) 
and for the CMB we use the 4-year COBE 53GHz map 
(full sky coverage.) The 31 GHz and 90 GHz companion 
maps both have significantly higher noise, as do the com- 
bination maps constructed by the COBE team to min- 
imize Galactic emission. Above a Galactic latitude of 
\bgai\ > 10°, the NVSS source counts are essentially all 
extra-Galactic. Prior to cross-correlating the two data 
sets, a secant law Galaxy model and a dipole (the largest 
structure in the 53 GHz map) were fit and removed. The 
effect of these corrections was minimal. 

The NRAO VLA Sky Survey (NVSS) is a 2.5 x lO'^ 
Jansky (2.5mJy) flux limited survey at a frequency of 
1.4 GHz It is complete for declinations 6 > -40° 

and contains > 1.8 x 10^ sources, with a mean source 
number density of 51.8 per square degree. While the dis- 
tances of individual sources are largely unknown, typical 
luminosity function models (e.g., |T5|Jl^ ]) indicate they 
are distributed in the redshift range < z < 2 with a 
peak in the distribution at z 0.8. This distribution, 
combined with the nearly full-sky coverage, renders the 
NVSS survey an ideal matter density probe with which 
to investigate the ISW effect. 

The NVSS survey was converted into a map using 
an equatorial quadrilaterized spherical cube projection 
[ p7t , the standard format for the COBE data. Those 
1.3° X 1.3° pixels that were only partially contained in 
the survey region were omitted. In making these maps, 
we discovered that the surface density of sources var- 
ied by ^ ±5% within several declinations bands. The 
more prominent of these bands coincide with discontin- 
uous changes in the rms noise levels in the NVSS survey 
. To correct for this, random sources were added to or 
subtracted from each pixel to eliminate the band struc- 
ture. The resulting map shows no declination dependent 
structure at level of < 1%. For comparison, the Poisson 
noise per pixel is ~ 11%. 

In order to exclude Galactic sources as well as nearby 
clusters of galaxies, the region within ±10° of the Galac- 
tic plane was removed as were 29 regions that con- 
tained pixels with source counts greater than 4cr (43%) 
above the average. This procedure cleaned the map of 
10 "objects" located more than 10° from the Galactic 
plane. Among these are the Orion Nebula and the nearby 
Virgo, Perseus, and Fornax clusters. That the moder- 
ately nearby, rich Coma cluster of galaxies was not one 
of the regions cut indicates to us that this windowing re- 
moved only Galactic and nearby extra-galactic sources. 

After the cleaning and correction operations, the map 
was repixelized by combining groups of four pixels into 
larger 2.6° x 2.6° pixels. The angular resolution of the 
COBE map is ~ 7° and this coarser pixelization is the 
standard format provided for the COBE DMR data sets. 
Each of the larger pixels was assigned a weight Wj cor- 
responding to the number of subpixels contributing to 



it and its value was multiplied by 4/wi to achieve the 
proper normalization. The radio source auto-correlation 
function (ACE) was computed according to 

i^BMiO) =Y,MN^ - N)^iiN, - N) /^WiWj, (1) 

where the sums are over all pairs of pixels with angu- 
lar separation 9, N is the mean number of sources and 
Wi is the pixel weighting factor. The measured ACE is 
displayed in Figure 1 with error bars that were deduced 
from Monte Carlo simulations. 
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FIG. 1. The auto-correlation function of the NVSS and 
the predictions from a ACDM theory, where Qa = 0.7 and 
h — 0.7. The inferred linear bias factor is bn = 1.6. The 
theory has been convolved with the pixel window function. 
Note that the errors are highly correlated. 

The effect of the above mentioned declination band 
corrections on the ACE was considerable. For 9 > the 
ACE of the corrected map is ~ 2-3 times smaller than 
that of the uncorrected map. This is not unexpected 
given the prominent band structure in the uncorrected 
map. To test the sensitivity of the ACE to our cor- 
rections, more aggressive band removal was effected by 
adding (and subtracting) sources at random from even 
more narrow declinations bands (8°, 4°, and 2°) in order 
to force the mean number densities in these bands to be 
equal. In all such cases, the ACE's were not significantly 
different from that of Figure 1. We conclude that our 
ACE is reasonably robust. 

One can determine the radio source bias factor by com- 
paring the measured ACE to the mass density correla- 
tion predicted theoretically assuming the COBE normal- 
ization and a given cosmological model. The curve in 
Figure 1 is the prediction of a sample ACDM cosmol- 
ogy convolved with window functions for the quad-cubed 
pixelization. The inferred radio bias parameter depends 
on the exact cosmology, as is shown in Figure 3, but we 
find typical values in the range 6^ = 1.3 — 1.6. Since we 
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compare on very large angular scales where the fluctu- 
ations are linear, the inferred biases are independent of 
the vagaries of non-linear structure formation. 

An essential ingredient in interpreting the auto- and 
cross-correlations is the source number-redsliift distribu- 
tion, dN/dz. The distribution used in Figure 1 was de- 
rived from a radio source luminosity function (LF) of 
Dunlop and Peacock (mean-z model 1) This LF is 
consistent with flux limited number counts of the present 
NVSS survey as well as with recent, deep redshift surveys 
p8[ . The mean-z model 3 of Dunlop and Peacock and 
the model of Condon are also reasonably consistent 
with the data and were also considered. The bias fac- 
tors deduced from these latter two LFs are within 15% 
of that deduced above, which indicate that the results are 
not overly sensitive to the details of dN/dz or the LF. 

We can compare our results to the ACF of the FIRST 
radio survey, which was determined by Cress et al. ||l9| ] 
and Magliocchetti et al. (20| toheujRR{9) ~ l-2xl0"3/6' 
for 0.01° < 6* < 4° where 9 is in degrees. While the 
FIRST survey is somewhat deeper (~ ImJy) than the 
NVSS survey, the luminosity function is such that the 
redshift distribution of these sources is nearly identical to 
that of the NVSS sources |lj]. Therefore, one expects the 
ACFs for these two surveys to be nearly the same. When 
converted to the the pixelization used in this paper, these 
results are consistent with that in Figure 1. 

The cross-correlation function (CCF) of the radio and 
CMB maps was computed using 
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iV)(r,-T)/af /^Wi/a^ (2) 



where Tj and cr| are the temperature and the noise vari- 
ance of the j*'* pixel in the CMB map, T is the mean 
CMB temperature, and the radio variables are defined 
as before. The COBE instrument noise is significant and 
the pixels were weighted accordingly in Eq. 2. However, 
if the pixels are weighted uniformly, the resulting CCF 
does not change significantly. 

While there may be unaccounted for systematics, it 
seems unlikely that these will be correlated in two such 
disparate maps. This is born out by the insignificant 
differences in the cross-correlation functions computed 
both with and without correcting the NVSS map for the 
declination band structure and with and without correct- 
ing the CMB map for both a secant law Galaxy and a 
quadrupole moment. 

The CCF is plotted in Figure 2 and is consistent with 
the data sets being uncorrelated. The effect of remov- 
ing the dipole and secant law Galaxy model from the 
CMB map was investigated by cross-correlating a set of 
Monte Carlo simulations in which a small correlated com- 
ponent having the same profile as the ACDM CCF was 
added to larger components that match the ACFs of the 
two maps. The results indicate that the dipole/secant 
law corrections significantly attenuate the CCF but only 



for separations, 6 > 15°. Therefore, these data were 
ignored in the subsequent analysis. The error bars in 
Figure 2 were computed from the Monte Carlo trials and 
are highly correlated. 
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FIG. 2. The cross-correlation of the NVSS with the COBE 
microwave background 53GHz map. The errors are highly 
correlated. 

The theoretical predictions for the cross correlations 
are calculated following the formalism of Refs. @,|l^. In 
order to simplify comparison to the data, we consider 
a single family of possible correlation functions having 
the shape of that in Figure 2 but with variable ampli- 
tude. This family closely matches the correlations pre- 
dicted for the range of ACDM cosmological models we 
consider. The 95% CL (1.65 a) upper limit for the zero 
lag correlation is w_r.t(0) < 0.038 while wrt{0) < 0.048 
at 98% CL (2.05 a.) 

The cosmological implications of this bound are sum- 
marized in Figure 3, where we plot the allowed region in 
the Hq — plane, where Qm is the fraction of the criti- 
cal density in both baryons and dark matter. Also shown 
on the plot are contours of constant radio bias and con- 
tours of constant F = flmh [h = Ho/100kms~^Mpc~^), 
which determines the shape of the matter power spec- 
trum. To first order, the predicted cross-correlation de- 
pends only on the value of the cosmological constant, 
and models with ^ 0.7 are excluded. However, 
there is some dependence on the Hubble parameter be- 
cause the radio auto-correlation includes contributions 
from smaller scales than generally contribute to the cross- 
correlation. Thus, models with more small scale power 
(large F, Hq) tend to have smaller biases and predict less 
cross-correlation than models with lower F and Hq- 

Our results have thus far assumed a constant linear 
bias factor, but in principle bias could be time and even 
scale dependent. Simple models of linear bias evolution 
lEll indicate that it is tied to the evolution of the lin- 
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ear growth factor D(t) (normalized to be unity today): 
bo — 1 = {b{t) — l)D{t). Thus, models that are presently 
positively biased were even more biased in the past. We 
do not expect this time evolution to change our results 
significantly, however. First, it is the levelling off of the 
growth factor D{t) at z ~ 1 that results in the ISW ef- 
fect. Thus the bias is also relatively constant during this 
epoch. Second, the effect of bias evolution is the same 
as changing dN/dz, and our signal is not very sensitive 
to such changes. Scale dependent biasing, if arbitrary, 
makes any predictions for large scale structure problem- 
atic. However, assuming the scale dependence is rela- 
tively weak, then so long as we are measuring the radio 
bias on roughly the same scale as the expected correla- 
tions, our results should be fairly robust. 
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FIG. 3. The region excluded at the 95% CL in the Hubble 
constant-matter density plane. Also shown are contours of 
constant P. The radio biases are inferred from the COBE 
normalization and the radio auto-correlation function. The 
boxed region represents the observationally preferred region. 

The integrated Sachs- Wolfe effect, if detected, would 
provide an important confirmation of cosmological the- 
ory and provide a new mechanism to distinguish between 
models with a cosmological constant from quintessence 
dominated or open universes. Of course, if radio sources 
themselves contribute significantly to the CMB then such 
a correlation could be explained otherwise. However, cur- 
rent models indicate that this is not the case. From flux 
limited 1.4 GHz source counts it can be shown that the 
contribution from these sources is less than 10% of the 
expected ISW contribution to ujut (0) and is even less for 
6* > 0. 

While the upper limit derived here is consistent with 
the currently favored Q\ ~ 0.65 CDM universe, it is only 
marginally so. It is worth noting that a similar limit has 
been derived from the frequency of strong gravitational 
lensing |22]. Since these observations are so close to be- 



ing in conflict with the leading cosmological model, it is 
important to pursue both types of investigations further. 

The cross-correlation constraint will be greatly 
strengthened by NASA's recently launched Microwave 
Anisotropy Probe (MAP), which will provide a full sky 
map with much lower noise and higher angular resolution 
than the COBE map. Monte Carlos simulations indicate 
that the cross-correlation of MAP with the NVSS survey 
should be able to detect the ISW effect for an Q\ = 0.6 
universe at the 95% confidence level. If, in addition, a fu- 
ture radio survey is able to increase the number of sources 
by a factor of ^ 3, then the Poisson counting noise will 
be effectively eliminated and the detection should be at 
the 3(7 level. Finally, if in the future there becomes avail- 
able a full-sky quasar or distant galaxy redshift survey 
with several million objects, one will be able to construct 
the 'ideal' redshift distribution function and approach the 
optimal 5.5cr detection derived by Crittenden and Turok 
for a Oa = 0.6 universe 0. These will constitute cru- 
cial tests which the ACDM universe must pass if it is to 
remain the favored cosmological model. 
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